Multi-scale fluctuations near a Kondo Breakdown Quantum Critical Point 
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We study the Kondo-Heisenberg model using a fermionic representation for the localized spins. 
The mean-field phase diagram exhibits a zero temperature quantum critical point separating a spin 
liquid phase where the f-conduction hybridization vanishes, and a Kondo phase where it does not. 
Two solutions can be stabilized in the Kondo phase, namely a uniform hybridization when the band 
masses of the conduction electrons and the f spinons have the same sign, and a modulated one 
when they have opposite sign. For the uniform case, we show that above a very small Fermi liquid 
temperature scale (~ 1 niK), the critical fluctuations associated with the vanishing hybridization 
have dynamical exponent z — 3, giving rise to a specific heat coefficient that diverges logarithmically 
in temperature, as well as a conduction electron inverse lifetime that has a T log T behavior. Because 
the f spinons do not carry current, but act as an effective bath for the relaxation of the current carried 
by the conduction electrons, the latter result also gives rise to a T log T behavior in the resistivity. 
This behavior is consistent with observations in a number of heavy fermion metals. 

PACS numbers: 71.27.+a, 72.15.Qm, 75.20.Hr, 75. 30. Mb 



I. INTRODUCTION 

A large number of experiments have been performed 
on metallic heavy fermion compounds close to a zero 
temperature phase transition (a quantum critical point 
(QCP)) driven by applied magnetic field, chemical dop- 
ing or pressure^"-^ . In the quantum critical regime, 
the thermodynamics and transport properties indicate a 
breakdown of the Fermi liquid. In many cases, the resis- 
tivity is quasi-linear in temperature over several decades, 
and the specific heat coefficient diverges logarithmically. 
The spin susceptibility typically exhibits an anomalous 
exponent in temperature. Neutron scattering experi- 
ments on some of these materials have revealed that the 
anomalous exponent in the dynamical susceptibility is 
identical for all points in the Brillouin zone^:^, suggesting 
a local character for the fluctuations. deHaas-vanAlphen 
experiments also find a divergence of the effective mass 
when approaching the critical point, along with a change 
in the Fermi surface topology when going through it^. 

These unusual observations have motivated many the- 
oretical studies that have attempted to capture these ef- 
fects. Most theories^i^'i'^i^ are based on the assumption 
that at the QCP, a spin density wave forms, and there- 
fore the critical fluctuations that destabilize the Fermi 
liquid are magnetic in nature^i^iii^i^. In three dimen- 
sions, these theories fail to capture simultaneously the 
linear temperature dependence of the resistivity, the log- 
arithmic divergence of the speciflc heat coefficienbiS,, and 
the anomalous exponent of the spin susceptibilitj^. For 
an antifcrromagnctic spin density wave transition, a cen- 
tral problem is that the critical fluctuations are confined 
to an inverse coherence length about the spin density or- 
dering vector, and consequently, only parts of the Fermi 
surface couple effectively with the critical bosonic modes. 

More recently, the problem has been approached from 



another perspective which takes the point of view that at 
the QCP, magnetic fluctuations suppress the formation 
of the heavy Fermi liquid, driving the effective Kondo 
temperature of the lattice (Tr-) to zer o^°i^^i^'^i^'^ . In 
this picture, the QCP is a fractionalized critical point 
at which the heavy quasiparticle deconflnes into a spinon 
and holon. One feature that distinguishes between these 
two classes of theories is that the flrst predicts the Fermi 
surface to change smoothly across the QCP, while the sec- 
ond predicts an abrupt change^-. Recent results of the 
Hall effect for YbRh2Si2i^, as well as the earlier men- 
tioned dHvA datai^, have lent support to theories of the 
second type. 

Here, we explore the possibility that in the quantum 
critical regime, the magnetic fluctuations are not the 
dominant ones at the QCP, and that the unusual behav- 
ior in thermodynamics and transport is due to critical 
fluctuations of a non-magnetic order parameter associ- 
ated with the vanishing energy scale Tk- One motivation 
for this point of view is the fact that in some compounds 
like YbRh2Si2, the gain in entropy inside the magneti- 
cally ordered phase represents only a few percent of the 
total entropy per localized spin^^. The order parameter 
we advocate is the field a associated with the hybridiza- 
tion between the localized spins and the conduction clcc- 
tronsilii^. At the QCP, the effective Kondo temperature 
for the lattice goes to zero, leading to a 'Kondo break- 
down' of the heavy Fermi liquid. The critical fluctuations 
of (J are gapless excitations, and we study how these fluc- 
tuations influence the properties of the metal using the 
formalism of the large N Kondo-Heisenberg model. 

There have been several earlier studies of this 
modeli^iii^ii^. Beyond the mean-field level, the Kondo- 
Heisenberg model can be treated as a lattice gauge the- 
ory. Senthil et. alM have examined the effect of the 
gauge fluctuations in this model, while Coleman et. alr^ 



studied the zero temperature transport anomalies. In our 
work, we find a number of novel effects associated with 
the fluctuations of the a field which were not discovered 
in these earlier studies. 

At the Kondo breakdown QCP, the metal passes from 
a magnetic phase (which we approximate, as in earlier 
work^'^, as a uniform spin liquid) to a Kondo phase. In 
the spin liquid phase, the f spinous are characterized by a 
'Fermi surface' which generically differs in size from the 
conduction electron Fermi surface. In the Kondo phase, 
these two surfaces become coupled due to the non-zero 
expectation value of a. In our study, we observe two new 
phenomena associated with this. First, for the case where 
the spinon and conduction electron masses have opposite 
sign, (T can order at a finite wavevector, leading to spa- 
tial modulations of the Kondo hybridization analogous 
to the LOFF state of superconductivit y^^i^^ . Second, we 
find the presence of multiple energy scales, spread over a 
very large range in energy, due to the mismatch between 
the two Fermi surfaces. The lowest scale, below which 
Fermi liquid behavior is restored, is extremely small (of 
order 1 niK), above which, up to an ultraviolet cutoff 
of order the single ion Kondo temperature, the critical 
fluctuations of a exhibit a dynamical exponent z = 3. 
This gives rise to a marginal Fermi liquid like behavior 
in d = 3 for the conduction electrons along the entire 
Fermi surface, due to scattering with the critical fluctu- 
ations. This property is to be contrasted with antiferro- 
magnetic spin density wave models, where only on parts 
of the Fermi surface the scattering of the electrons with 
the critical mode is effective. Next, since the f spinous 
do not carry current, but act as an effective bath for the 
relaxation of the current carried by the conduction elec- 
trons, the marginal Fermi liquid behavior also gives rise 
to a resistivity that goes as TlogT. This behavior is un- 
like either that of ferromagnetic spin density wave models 
in which the transport lifetime is less singular than the 
single particle lifetime (i.e., in the latter models, forward 
scattering does not degrade the current), or that of an- 
tiferromagnetic spin density wave models in which the 
"cold" parts of the Fermi surface dominate the trans- 
port propertiea^^. Moreover, a logarithmic dependence 
is found for the speciflc heat coefficient from both the 
gaugeiS, and a fluctuations. The latter also give rise to 
an anomalous temperature exponent of 4/3 in the uni- 
form spin susceptibility. A summary of our results have 
been presented in a shorter paper~ . 

The phenomenon of the breakdown of the Kondo effect 
at a QCP can also be studied in the more general con- 
text of a periodic Anderson model. This generalization 
is discussed in other works2ii2^. 



II. MODEL AND FORMALISM 

The starting point of our theory is the microscopic 
Kondo-Heisenberg model in three dimensions, which de- 
scribes a broad band of conduction electrons interacting 



with a periodic array of localized spins through antiferro- 
magnetic Kondo coupling Jk > 0. Additionally, the lo- 
calized spins interact with one another via nearest neigh- 
bour exchange Jh > 0. The Hamiltonian for the large 
N version of this model, where N denotes the enlarged 
spin symmetry group SU{N), is given by 



{ij),a i,a,P 

Jh 



+ AT 2^ JiaJit3Jjf3JJa- 



(1) 



{ij),a,P 



Here c\^ {cia) are creation (annihilation) operators for 
the conduction electrons with spin index a — {1,N) at 
site i, and (ij) refers to nearest neighbour sites, t is 
the hopping matrix element between neighbouring sites 
for the conduction electrons. The SU{N) generalization 
of the locahzed spins 5° with a = (1, ..., 7V^ — 1) at each 
site i are expressed in terms of Abrikosov pseudofermions 
(or spinous) by Sf = E„/3 fLiKf,/^)U, where F" are 
the generators of the SU{N) group in the fundamental 
representation. This fermionic representation of the spin 
operator gives rise to a local constraint at each site, i 



Tf' f = - 



Wi. 



(2) 



We note that in the context of the heavy fermion sys- 
tems, the Heisenberg exchange term is often equated to 
the RKKY interaction between the localized spins which 
is mediated by the mobile conduction electrons. In such 
a scenario, the Heisenberg coupling Jh (x. PqJ]^, where 
po is the density of states of the conduction electrons at 
the Fermi level. However, for the purpose of the present 
study, it is convenient to consider Jh as a parameter in- 
dependent of Jk- Microscopically this can be justifled by 
noting that, in principle, there can be other sources which 
generate this coupling, such as superexchange within the 
narrow band of /-electrons. 

In order to perform a systematic large TV study of 
the system defined by Eqs. ^ and ([2]), the first step 
is to decouple the interaction terms which are quartic in 
fermionic operators using a Hubbard-Stratonovich trans- 
formation. The Heisenberg exchange term is decoupled 
using a bosonic link variable 0y — > ^^ fJafja, while the 
Kondo interaction is decoupled by introducing a complex 



bosonic field a, 



y f 



In the next step, follow- 



ing Ref. 13, we assume that in three dimensions, 0y con- 
denses in a uniform spin liquid phase, i.e., {4>ij) — 4>q at 
the mean field level. This provides a dispersion to the 
spinon band which, as we show later, is an essential in- 
gredient to obtain the breakdown of the Kondo effect. 
We note that there is no clear evidence of a spin liquid 
phase in any heavy fermion system near its quantum crit- 
ical point. Rather, the typical phase diagram exhibits a 
QCP that separates a magnetic ground state (typically 
an antiferromagnet) from a paramagnetic heavy Fermi 



liquid. Consequently, it is useful to discuss the motiva- 
tion for our choice of a uniform spin liquid phase for the 
Heisenberg link variable (pij . This choice is partly guided 
by convenience: since our main purpose is to study the 
consequences of the breakdown of the Kondo effect, the 
choice of a uniform spin liquid can be viewed as the sim- 
plest device which allows the vanishing of the Kondo en- 
ergy scale (indicating the breakdown of the Kondo effect) 
at the mean field level. More physically, one can view 
the uniform spin liquid as a mean field description of 
the short range magnetic correlations that persist when 
a magnetic ground state is destroyed by quantum fluctu- 
ations. However, to demonstrate this point concretely is 
not simple, and beyond the scope of the present study. 
The key point is that the spin liquid provides a band- 
width for the f electrons. Other approaches, for instance 
one where the bandwidth is due to direct f-f hopping, 
should yield similar results in regards to the breakdown 
of the Kondo effect that we describe here. 

The system can now be described by the Lagrangian 



^ = E [4(5.%-i)c,a + /l((a.-A,)5, 
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0oe"")/j"] + 



^V-^ ^ V- t 



Jh 



E 



y^iaJia'^i 



H.c. 



(3) 



where V (the volume of the system) is set to 1. In the 
above, Xi are Lagrange multipliers (scalar potential) that 
enforce the local constraint of N/2 spinous per site. Now, 
given a many-body wave function that satisfies this con- 
straint, a single hop of a spinon takes the state out of 
the physical subspace. Consequently, for the kinemat- 
ics of the spinous, only simultaneous opposite hops be- 
tween two neighbouring sites is a physically allowed pro- 
cess. This implies that the local spinon current operator 
Jfi = at every site i. The gauge fields a^ (vector po- 
tential), associated with the phase of 4>ij , ensure that this 
constraint is satisfied. The appearance of the scalar and 
vector potentials can also be understood by noting that 
C is invariant (up to a term which is a total derivative 
of imaginary time) under a local U{1) gauge transfor- 
mation /iQ, -> /laC'^S CTi -^ CTjC"*^', Xi ^ Xt + idrOi, 

ttij -^ Qij — 9i + 9j, a consequence of the fermionic rep- 
resentation of the spin and the constraint Eq. ([2|— . 

In the following we examine the above Lagrangian, first 
in a mean field approximation and then consider Gaus- 
sian fluctuations of the action around the mean held solu- 
tion. This involves studying the possibility of hybridiza- 
tion between the conduction and the spinon bands (for 
(ai) ^ 0) as well as calculating the hybridization fluctu- 
ation which is an interband particle-hole excitation. As 
such, one needs to characterize the dispersions of the con- 
duction and the spinon bands. We do this by assuming 
that the bands have a parabolic dispersion (to facilitate 
calculations) , and we introduce the following two impor- 
tant parameters. First, a = 0o/-D, is the ratio of the 
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FIG. 1: (Color online) Dispersion of conduction and spinon 
bands, with the mismatch wavevector, q* , and the mismatch 
energy. Ex = avpq* , indicated, where a is the ratio of the 
spinon and conduction bandwidths. An artificially large value 
of a was used in this plot (0.5) so as to better illustrate the 
origin of E^- 



spinon bandwidth ^q and the conduction bandwidth D. 
As we will see in the next section, at the Kondo break- 
down QCP 00 ~ Jh ~ ^0., where T°. = De'^/'^P"^'^'' is 
the single-ion Kondo energy scale of the system, which 
is typically of order 10 K in heavy fermion systems. As- 
suming D ^ 10'* K, we get a ~ 10"'^. Second, while 
the spinon band is half filled due to the constraint (for 
A^ = 2), the conduction band filling is generic. With- 
out any loss of generality, we take the conduction band 
to be less than half filled. This implies that the Fermi 
wave vector of the conduction band kp is different from 
that of the spinon band fc^o- We denote this mismatch by 
q* = kpQ — kp, and assume that the fraction {q* /kp) is of 
the order 0.1—. This would mean that while kp and kpo 
are of the order of the Brillouin zone dimension, the mis- 
match wave vector q* is one order of magnitude smaller. 
The parameters a and {q* /kp) affect the important en- 
ergy scales of the system. This is illustrated in Fig. [1] 
where we show the conduction and spinon dispersions. 



III. MEAN FIELD TREATMENT 

At the level of the mean field approximation, we re- 
place the bosonic Hubbard-Stratonovich fields and the 
Lagrange multipliers by their expectation values, and we 
study the approximate Lagrangian given by 

Cmp = E H"(^^'^y ~*)^J" + '^»"((^^"^-^*))'^*J 

{ij)a 

- '^o)/ja] + YE<^^) + j^El<'^«)l +^7^ 

i i 

+ E(4/»(^^) + H.c.). (4) 



In the following we write the dispersion (et) of the con- 
duction band as 



.k = .+ -^, 



(5) 



where e = vpik — kp), k is the magnitude of k, and vp 
is the Fermi velocity of the conduction electrons. The 
dispersion {e'^) of the spinon band is similarly written as 



4 



a 



{e^VFq*) + 



*\2 



D 



(6) 



We note that, in the above, both the bands are taken as 
electron-like, for which we find that the mean field equa- 
tions yield a spatially uniform solution, namely, {ai) = Uf) 
and (Ai) = Aq. In the case where one of the bands is cho- 
sen to be hole-like, we find a spatially modulated solution 
which we discuss in appendix I A II The free energy cor- 
responding to Eq. ^ is given by 



Fmf 

N 



a 



Tr [ln(-G-i(iL^„,k)) +ln(-G,-i(ic^„,k))] 



'^l , ^0 



Jk Jh 



(7) 



where (3 is the inverse temperature, ujn is the fermionic 
Matsubara frequency, and Tr corresponds to a trace over 
space-time co-ordinates. In the above 



•^a.U*^"'^) = *^" 



a.h 



where 



a.b 



Ek 
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^ir 



4a2 



(8) 



(9) 



We evaluate the free energy given by Eq. ([7]) at zero tem- 
perature (T = 0) in the limit {q* /kp) -^ 0. The details 
of this evaluation is given in appendix I A 21 As a function 
of a and co: f^nd to 0{a'^) accuracy, we find {a <^ 1) 



F, 



MF 



2 r 



N 



PoD 



ln(l 
a 



2poJh 






Poo-Q 



const, 



PqJk 



(10) 



where the constant part has explicit Aq dependence. 
Since the precise value of Aq is of no importance for our 
results, in the following we ignore the mean field equa- 
tion for Aq. Minimizing Fmf with respect to a and ctq 
we get 



0, (11) 



PoD^ 


( a 


1\ 2ol Aat 1 
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\poJh 


2poao 


\f ^ 


A^y^fo] 


\poJk 



= 0, 



(12) 



respectively. We study these equations by keeping the 
Heisenberg parameter Ju fixed, while varying the Kondo 
parameter Jx, and find two solutions corresponding to 
two mean field ground states, (i) First, a uniform spin 
liquid phase where ctq = 0, which implies that in this 
phase, the Kondo effect fails to occur and the localized 
spins remain unscreened in a uniform spin liquid state. 
In this phase, a = ag = (po>^if)/3, which implies that 
the Heisenberg coupling sets the scale for the spinon dis- 
persion, since 4>q = {pqDJh)/6 '^ Jh- It is simple to 
check that this solution is stable for Jk < Jk^ , where 



ln(- 
PqJk. \ a 



(13) 



(ii) For Jk > Jk^ the stable mean field solution corre- 
sponds to ao ^ 0, indicating a ground state where the 
local moments are screened by the Kondo effect and a 
heavy Fermi liquid is established below an energy scale 
Tk ~ TrpoCQ. The growth of the Kondo order parameter 
in this phase is given by 



CTo oc Jh In 



Jk — Jxa 



1/3 



D 



(14) 



where /3 = 1/2 is the typical mean field exponent. We 
also find that the spin liquid order parameter decreases 
in this phase, and is given by 



= ao-^-^ + 0{4). 



(15) 



Thus, from the above mean field study, we find that, 
in the presence of a finite bandwidth of the spinous, the 
Kondo effect takes place only when the Kondo coupling 
Jk is larger than a finite value Jk^ ■ This establishes the 
Kondo breakdown QCP where the lattice Kondo energy 
scale Tk vanishes. In the current formulation of the mean 
field theory, the Kondo breakdown QCP separates a uni- 
form spin liquid ground state {Jk < Jk^) from a heavy 
Fermi liquid ground state {Jk > Jk^)- It is important 
to note that if we define a single-ion Kondo scale (T]^-) 
as a function of Jk for the system by 



T%{Jk) = De 



-1/(P0./k) 



(16) 



using Eq. P^ we conclude that at the QCP 

Jh-T^k{Jk:). (17) 

This shows that the Kondo breakdown QCP is estab- 
lished as a result of a competition between the Kondo 
energy scale and the magnetic energy scale, even though 
there is no long range magnetic order in the present 
study. The reduction of the spin liquid order parame- 
ter, given by Eq. p^ . provides further evidence for this 
competition. Therefore, this mean field study can be 
viewed as a microscopic realization of the energetic argu- 
ment that Doniach had proposed several decades ago for 
the existence of a QCP in heavy fermion systems^^. 



IV. FLUCTUATIONS 

In this section, we study the niassless fluctuations in 
the quantum critical regime. There are two such modes: 
(a) one associated with the phase of (pij which are the 
gauge fluctuations, and (b) the fluctuations of the com- 
plex order parameter (cr|,cri) which are gapless due to 
the vanishing of the Kondo energy scale Tr- at the Kondo 
breakdown QCP. 



A. Gauge Fluctuations 

Since the gauge fluctuations of the system have been 
studied earlier—, here we just summarize the main 
points for the sake of completeness. It is convenient 
to work in the Coulomb gauge V • a = 0, where the 
vector gauge fields a^ (/i = x, y, z) are purely trans- 
verse^. In this gauge the fluctuations of the scalar po- 
tential A decouple from a^, and give rise to a screened 
Coulomb interaction between the spinous which can be 
neglected. Next, since the fields a^ enter the the- 
ory as vectorial Lagrange multipliers to satisfy the con- 
straint that the local spinon current is zero, they be- 
have as 'artificial photons' without any intrinsic dy- 
namics of their own. Their dynamics is entirely gen- 
erated by their coupling to the matter field, namely 
the spinon band, and therefore these bosonic modes are 
overdamped. The propagator for the transverse gauge 
fields is defined as £)^i,(x,r) = (T.^ [a^(x, T)ai/(0, 0)]), 
which in frequency-momentum space has the standard 
form D^,y{c\,iVln) = ((5^^ - qf_,qu/q^)'n~^{q,i^n), with 
n((7, ifi„) oc [{q/2kpQ)'^ -j- \Q,n\/{avpq)]. Here n„ is a 
bosonic Matsubara frequency, and the above expression 
for the gauge propagator D^i,(q, ifi„) is valid for fre- 
quencies smaller than the spinon bandwidth aD. As a 
result, the gauge excitations are characterized by a dy- 
namical exponent z — i, which in rf = 3 are known2^ 
to give a contribution to the specific heat coefficient 
7 = —d'^F/dT'^ oc \n{aD/T) and to the static spin sus- 
ceptibility 5xs fx T^h\{aD/T). Finally, it has been ar- 
gued in the literature that the gauge fluctuations convert 
the flnite temperature mean field phase transition line 
into a crossover linei^i^. 



B. Fluctuations of the Kondo Boson 

At the QCP, where the Kondo coupling is tuned to 
its critical value J^^, the critical fluctuations of the 
continuous phase transition are given by those of the 
complex order parameter fields {a'^ , a) . The propa- 
gator for these fluctuations is defined by D„{x,t) = 
{Tr [at(x,T)a(0,0)]). We get i^-i(q,ifln) = 1/Jk + 
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FIG. 2: (a) Vertex for the interaction between the con- 
duction electrons (solid line) and the spinons (dashed line) 
mediated by the hybridization fluctuations a (wiggly line), 
(b) The interband polarization involving conduction electrons 
and spinons, which generates the dynamics of a. For momen- 
tum transfer q > q* , where q* is the mismatch between the 
conduction and the spinon Fermi surfaces, a is an overdamped 
critical mode with dynamical exponent z = Z. 



n/c(q,«il„), where 



n/c(q, ifi„) 



E 



Gc(k, iw„)G/(k-}-q, iw„ - iUn) 

(18) 
is the interband polarization bubble between the conduc- 
tion and the spinon bands. In the above G'^^{}i.,iujn) = 
{iojn — Ck) is the propagator for the conduction elec- 
trons, while G'7^(k, iw„) = (iw„ — e^) is the propaga- 
tor for the dispersive spinons. We write n/c(q, if^n) = 
n/c(q, 0) + An/c(q, if^n), where n/c(q, 0) is the static 
part of the fluctuations and An/c(q, i^n) is the dynamic 
part. We first compute the static part which can be writ- 
ten as 



n/c(q,o) = 5] 



nF(ek)-nF(e£+q) 



Ek - e 



(19) 



k+q 



where npi^) is the Fermi function. We find that n/c(q, 0) 
is independent of momentum if the dispersions are lin- 
earized in Eq. ()19p . This implies that the momentum de- 
pendence is due to k r^ kp in the k-integral of Eq. (|19p . 
for which it is important to retain the quadratic disper- 
sions of the bands. Furthermore, since the main con- 
tribution is for fc ~ fcf, the small momentum scale q* 
is unimportant and can be set to zero to facilitate the 
calculation, and we write ek = (fc^ — kp)/(2m) and 
^k ^ (^"^ ^ fcF)/(2»Tio). Then, in terms of a, the ratio of 
the two bandwidths that we introduced earlier, we have 
a = m/mQ. Using k <-> k-l-q inside the k-summation we 



get 



where 



n/c(q,0) 
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C + qk/niQ 



where B — A{k'^ 
g^/(2mo), and A 



C — qk/mo 

- fc2) - gV(2™), C . 
[l/(2m) - l/(2mo)]. 



After perform- 



ing the momentum integration, we expand the result- 
ing expression in powers of {q/kp), and we use po = 
rnkp/i'^T^'^), the density of states per spin of the conduc- 
tion electrons at the Fermi level. To leading order in 
(q/kp) we get 



n/c(q,0) = po 



Pa 



In a 1 — a^ 


-|- 2a In a 


fq 


a -a) 4(1 
- li^ - + 77T 


-a)3 


ykp 



(20) 



Note that the po ln(a) term in the above equation has 
been derived in appendix A2 using a slightly different 
method for the calculation of the mean field free energy 
in Eq. (|10p . This term, along with 1/Jk, define the mass 
(1/Jif -|-Po1ii(q;)) of the Kondo boson, which goes to zero 
at the QCP. 

Next we calculate the dynamic part of the fluctuations 
which can be written as 

An/,(q,if]„) = n/c(q,iaO-n/c(q,o) 

= a y^ Gcjk, iuin) [Gf{k + q, iuin - i^n) 



G/(k-Kq, zcj„)] 



(21) 



Unlike in the case of the static part, here the dom- 
inant contribution is from the interband particle-hole 
excitations around the two Fermi surfaces, for which 
the spectra can be linearized. We write ek = e for 
the dispersion of the conduction electrons, and e^+q = 
a{e—VFq* +VFqz) for the dispersion of the spinous, where 
z is the cosine of the angle between wavevectors k and 
q. Approximating the k-summation by 



E 



Po 
2 



de I dz. 



at zero temperature we get 
Pa 



An/c(q, i^Ti 



2(1 -a) 



[Yi +Y2 + Y3 + Yi] , (22a) 
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El \ 

1 T ln(£'i T vpq) 

vpqj 



with El = vpq* — ifln/a, and 



Y 



3,4 



1 =F -^ ] ln(^2 T vpq) 
vpqj 



(22b) 



(22c) 



with E2 — Vpq* — ifin- From the above expression of the 
dynamic part given by Eqs. (|22ap - (|22cp . we next extract 
the leading behaviour in different regimes of frequency 
and momentum. For this we need to compare the mo- 
mentum q with q* , and the frequency f2„ (a continuous 
variable at T = 0) with the energy scales E^ = avpq* 
and avpq. Note that vpq* ~ 10^ K is an energy scale 
much larger than the ultraviolet cut-off of the theory 
aD ^ 10 K (the spinon bandwidth), and therefore we 
need to consider only |51„| <C vpq* . We find five distinct 
regimes which are as follows: 
(i) |f2„| < Ex and q < q* , where 



E 



Y 



-2(1 



Ex 



+ J(i + ") 



11 Lfi 

Ex 



(23a) 

Note that in the above, we retained two sub-leading terms 
because there are regimes where the sub-leading terms 
are larger than the static {q/kpY term, 
(ii) \Vln\ < Ex and q> q* , where 



^r,«-2(l-a) 



avpq 



i-sgn(0„) H 



(23b) 



(iii) |51„| > Ex and q < q* , where 
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E^^^ 



In 



Ex 



6 \q* 



Ex 



(23c) 



(iv) |ri„| > avpq > Ex and q > q* , where 



E^^ 



In I ^j +l + i-sgn{Q„) 



\avpqj 2 

(v) avpq > \fln\ > Ex and q> q* , where 

2i L'fi 



(23d) 



^r,«(l-a) 



avpq 



-i'KSgn{Qn) + {1 + a)- 






avpq 

(23e) 
At the quantum critical point, the mass of the Kondo 
boson goes to zero due to Eq. p^ . The leading fre- 
quency and momentum dependences of ^'^(q, «f^n) are 
determined using Eqs. ^, ^ and ^. The details of 
the various asymptotic structures of Z3o-(q, ^^n) in differ- 
ent regimes of frequency and momentum are discussed in 
appendix IB II Among the forms of D^lq^ifln) given in 



Eqs. (jBip - (|B3|) . only the following two asymptotic struc- 
tures are important for obtaining the leading contribu- 
tion of the Kondo boson to thermodynamic and transport 
properties. 

First, for |ri„| < [aD /{2TT)]{q* /kpY and q < q* , we 
get 



D-Hci,i^n) 



Po 



4 [kp 



2i L^i 



(24) 



which gives rise to an undamped propagating mode with 
dynamical exponent z = 2 (the dispersion of which is 
given by setting Eq. (|24p to zero) . The existence of this 
mode is a direct consequence of the mismatch between 
the Fermi surfaces of the conduction and the spinon 
bands. Due to this mismatch, a minimum momentum of 
q* is necessary to excite an interband particle-hole pair. 
Consequently, for momentum q < q* , the spectrum of the 
Kondo boson lies outside the continuum of the interband 
particle-hole excitations and thereby remains undamped. 
Note that this massless mode corresponds to hybridiza- 
tion fluctuations about the QCP, and becomes massive 
for Jfc < Jkc (this is realized by adding a constant term 
S to Eq. (pi)l ). Since Hfc at q=0 diverges logarithmically 
at i?a;, the mode energy never exceeds E^- The mode 
dispersion, which is quadratic about (7=0, is more com- 
plicated as q approaches q* due to logarithmic corrections 
to Ufa and is described in greater detail in appendix lB 21 
Second, for most of the phase space, the spectrum for 
the fluctuations of a lies within the interband particle- 
hole continuum, and we get 



D^\q,inn) ~ Po 



kp 



IT |ri„| 



2 avpq 



(25) 



i.e., an overdamped critical mode with dynamical expo- 
nent z = 3. Next we note that, since we assume q* <C kp, 
the overdamped z = 3 critical mode occupies most of the 
momentum space and therefore almost always it provides 
the leading contribution to thermodynamic and trans- 
port properties. In this regime, the scaling of frequency is 
given by Sl,i ~ [{aD) / {2'K)]{q / kpY , and since this regime 
ends for q < q*, one obtains the infrared energy scale 



E* 



caD{ ^ 

kp 



(26) 



where c is l/(27r). The true value of c is slightly smaller 
[^ 0.1) since there are logarithmic corrections to life as 
q approaches q* . A more detailed account is given in 
appendix lB 21 We note that E* , which can be estimated 
to be ~ 1 mK, appears as an infrared crossover scale for 
any physical property that is affected by the excitations 
of a. On the other hand, the ultraviolet cutoff scale is 
provided by aD ~ 10 K, which is the bandwidth of the 
spinous, or equivalently the single ion Kondo scale by 

Eq. inD. 



V. THERMODYNAMICS OF THE KONDO 
BOSON 

In this section, we study the effect of the fluctuations 
of the Kondo boson a on the thermodynamics of the sys- 
tem in the quantum critical regime. In particular, we 
compute (a) the contribution to the free energy, (b) the 
temperature dependence of the static spin susceptibility, 
and (c) the crossover lines in temperature which demar- 
cate the quantum critical regime. 



A. Free energy 

The contribution of the fluctuations of a to the free 
energy (per unit volume) is given by 

^^ f°° dVl / n\ r 1 

F = J2j^^—coth i^—j Imln [D-\q,n + i7j)] , 

(27) 
where Da{q, ^ + irf) denotes the retarded propagator for 
the Kondo bosons. We find that, for all temperatures 
T < aD, the leading T dependence of the free energy 
F is given by that part of phase space where the mode 
is overdamped (with dynamical exponent z = 3) and for 
which the expression for the propagator is approximately 
given by Eq. ([25]) . The details of this demonstration, 
as well as the evaluation of the sub-leading contribution 
from the other regimes, is given in appendix IB 31 For 
T > E* , the leading T dependence of F is given by 



dqq'^ 



(28) 



^-^^ / dficoth 

27r3 Jo 


v2r 


_i /27rn\ 




tan — ^ 

V r J 





Here q and qc are dimensionless momenta in units of kp, 
and f2 and T are dimensionless energies in units of aD. 
Since the g-integral is ultraviolet divergent, we use the 
Fermi momentum as an upper cutoff. The infrared cut- 
off, qc, for the z=3 regime is dependent on the particular 
temperature range considered, since the leading T de- 
pendence comes from frequencies i7 ^ T. For T > E* , 
q > rj^'^, for which we can approximate tan~"'^(a;) ~ x 
and replace the cutoff qc by Vl^/^ . Performing the inte- 
grals, we find 



F{T) 



1,3 



■■■(^ 



aD' 



T> E*. 



(29) 



We note that this contribution adds to a similar T^ ln(r) 
contribution from the transverse gauge fluctuations 
(which are massless z = 3 excitations). They give rise 
to a ln(r) behavior for the specific heat coefficient. 

For T < E* , the leading contribution to the free energy 
is again given by Eq. (^5)1 with qc = q* /kp for the infrared 
cutoff of the q- integral. This is because for Q, ^ T < E* , 
the z = 2) regime exists for q > q*. As a result, because 
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FIG. 3: (Color online) Numerical estimate of the contribution 
to the specific heat coefficient, C/T, coming from the Kondo 
boson. A value of a = O.OOf and q* /kp = O.f was assumed, 
with the momentum integral cut off at kp and the frequency 
integral at O.lEx- Note the logarithmic behavior of the z—3 
contribution which is cut off for T < E* (with E* ~ O.OOIB^), 
and the sub-leading nature of the z=2 contribution. 



J71/3 ^ q* /kp in this temperature regime, the lower cut- 
off remains at q* /kp. This gives. 



F{T) 



-4)K^ 






T <E* 



(30) 



This T^-dependence cannot be distinguished from ordi- 
nary Fermi hquid corrections, and in this temperature 
regime the free energy is dominated by the T^ ln(T) con- 
tribution from the transverse gauge fluctuations^^. 

The collective mode gives a magnon-like contribution 
to the free energy {F '^ T^/^), and is sub- leading rel- 
ative to the 2,-2) contribution (see appendix IB 3P . We 
illustrate this by showing in Fig. [3] a numerical deter- 
mination of the contribution of the specific heat coeffi- 
cient, C/T, coming from the Kondo boson, using the fc 
polarization bubble of Eq. (|22ap . In this plot, one sees 
the sub-leading contribution arising from the z=2 region, 
the logarithmic contribution from the z=3 region which 
saturates for T < E* , and the small difference between 
the positive and negative fl contributions from the z=3 
region due to the chirality of the fc polarization bubble. 



B. Static spin susceptibility 

At the mean field level, where the critical ffuctua- 
tions of a are ignored, the temperature dependence of 
the static spin susceptibility Xs{T) is entirely analytic, 
namely a constant (Pauli susceptibility) plus a T^ term, 
which is usual for band fermions. Next, when we take 
the critical fluctuations into account, we expect the cor- 
rection to Xs {T) to be non-singular (since the transition 



is non-magnetic and the excitations of a are in the sin- 
glet channel), but non-analytic (due to the massless ex- 
citations). In order to evaluate this temperature depen- 
dence, we first need to compute Da^(q,iCln) in the pres- 
ence of a magnetic field (B). For a finite B, the effect 
of the Zeeman term is to shift the Fermi wave vectors 
kp and kp of the conduction and the spinon bands, re- 
spectively. We get, kp^kpzt {fj,BgfB)/{avp), and 
kp ^ kp ± {^BgcB)/vp, where gf and gc are effec- 
tive Lande g-factors of the spinous and the conduction 
electrons, respectively, /is is the Bohr magneton, and 
± refers to the up and down spins, respectively. Since 
a <C 1, and in general gf > gc, we can ignore the cou- 
pling of B to the c-electrons and consider the effect of 
the Zeeman term as a renormalization of the mismatch 
wave vector q* , which is given by 



q*± 



iJ-BgfB 
avp 



Next, we note that, in the presence of a finite q* , one 
expects n/c(0,0) to have corrections of the type q* /kp 
and [q* /kp)^ (which are not calculated in Eq. (|20p . since 
the evaluation was performed in the limit q* ^0). This 
implies that, in the presence of a magnetic field, we ex- 
pect a correction to njc(0,0) which is proportional to 
[[^jLBgiB)/{aD)\'^ (since the excitation of a is in the 
singlet channel, we do not expect a linear term in B). 
Adding such a term to Da{(\,i^n), and noting that 
the leading temperature dependence is due to the over- 
damped z — i mode, we can generalize Eq. (|28p to obtain 
the B dependence of the free energy as 



F(S,T) 



k%aD 



X tan 



dri coth 


27rfi 



2r 



dqq 



h^q 



(31) 



Here energy and momenta are in dimensionless units (as 
in Eq. pSjl ) and h = {iJ,BgfB)/{aD) is the dimensionless 
magnetic field. Writing the correction to the static spin 
susceptibility due to the fluctuations of a as 6xs{T) = 
-[d^F/{dBf]B=o, we get for T > E* 



5xs{T)^-{^iBgf? 



24/3r(4/3)C(4/3) 

7^5/333/2 '^F 



rpi/3. 



{aDy/3 ' 
(32) 
while for T < E* , the lower cut-off is at q* , making the 
mode effectively massive, and we get 



Sxs{T) « -ifj-Bgff 



3{q 



1,5 
t\,p 



*\2 



{aDf 



(33) 



As in the case of the free energy, the non-analyticity in 
the leading temperature dependence is cutoff below E* 
due to the mismatch wave vector q* . As noted before, 
the gauge bosons give rise to a T^ ln(r) contribution to 

Xs- 



C. Crossover lines defining the quantum critical 
regime 

The crossover lines in temperature that demarcate the 
quantum critical regime are symmetric about the QCP 
S = Sc = 0, where S = ^/(poJk) — ^/(pqJkc) is the 
dimensionless tuning parameter of the theory (for fixed 
Jh)- On the heavy Fermi liquid side of the QCP, such 
a line usually defines the boundary of the finite temper- 
ature phase transition. However, it has been argued in 
the literature that the gauge fiuctuations convert the fi- 
nite temperature mean field phase transition line into a 
crossover lin o^^i^" . These lines are determined by the 
temperature dependent mass Sm{T) of the excitations 
of a. In a Ginzburg-Landau approach, these excitations 
are generated by the quartic uq IctI** coupling in the action, 
where uq ~ po/{ciD)^ from Eq. pU]) . In the following, we 
compute Sm{T) generated due to the propagating mode 
with z = 2 given by Eq. ([24]) , as well as that generated by 
the overdamped mode with z = 3 given by Eq. (|^ . The 
contributions from the other regimes of D^{q,iD,n) are 
always sub-leading. The general expression for 6m(T) is 
given by 

(34) 
Denoting the contribution of the z = 2 mode as 6mi{T), 
we get using Eq. 









where nB{x) — (e^^^ — 1)^^ is the Bose function. For 
the leading T-dependence, we write nh{x) ss l/x, with 
an appropriate ultraviolet cutoff for the g-integral. For 
T < E*, this cutoff is kpiT/E^y^^, and for T > E* , this 
cutoff remains at q* . We get 



Sra,m.(i^)lZ. T<E- 



f2uok'j,q* 



V T^^Pa 



(35) 



T, T> E* 



Next, denoting the contribution of the z = 3 mode as 
Sm^iT), we get 



5m2 (T) 



f AuoaDk^, 



\ tt'' 



00 



Po / JO 
„3 



dVl coth ( ^ ) ^ 



X 4^ '^^^FTi^' 



where q and q^ are dimensionless in units of kp^ and fi 
and T are dimensionless in units of aD. For T > E* , we 
can put qc '^ 17^^'^ for the leading term, while for T < E* 




FIG. 4: Self-energy diagram for the conduction electrons due 
to scattering from the critical excitations of a whose dynami- 
cal exponent is z — 3. In three dimensions, this has a marginal 
Fermi liquid form. 



we have qc 


= q* . This gives. 




5m2[T) « 


'2uok\,' 
. 3po . 


\q* ) aD' 




r27/3r(4/3)C(4/3)woA:?.1 


7-4/3 




33/27r5/3po 


(Qi?)l/3 



T> E*. 
(36) 



Comparing Eqs. ^ and §^ we find that, for T < E* , 
the leading T dependence is given by the z — 2 mode and 
Sm{T) fa SniiiT), while for T > E* , the leading term is 
from the z — 3 damped mode and 5m{T) « 5m2{T). 
Consequently, the crossover lines in temperature which 
define the quantum critical regime are given by 



Toe |<5 
o^\5 



,2/3 

I : 

,3/4 



T <E* 
T> E* 



(37) 



VI. 



QUASIPARTICLE LIFETIME AND 
TRANSPORT 



In this section, we first evaluate the quasiparticle life- 
time {tc) of the conduction electrons due to scattering 
from the excitations of a, and then argue that this life- 
time can be identified with the transport lifetime {rtr) 
for the evaluation of the temperature dependence of the 
resistivity. 

For the process shown in Fig. [51 the general expression 
for the imaginary part of the self-energy of the conduction 
electrons is given by 

/°° rlfl 
— [nsin) + npin - u;)] 
q -00 I" 

xlmD^iq,n)lmGfik-q,uj-n), (38) 

where np{x) = (e^'-^ -I- 1)~^ is the Fermi function, and 
R denotes retarded functions. At zero temperature, this 
gives 

ImEf(k,c.) = 5] r^Imi?^(q,f7) 
„ Jo I" 



xlmGf(k-q,u~n). 



(39) 
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We evaluate the above expression for a conduction elec- 
tron on the Fermi surface, i.e., for |k| = kp, and we find 
that the leading frequency dependence is always due to 
the overdamped z = 3 mode whose expression is given 
by Eq. ((25)) . The z ~ 2 mode does not contribute since it 
cannot kinematically connect the / and c electrons. We 
write 



E 



1 

4^ 



dqq'^ 



dz, 



where z is the cosine of the angle between k and q. After 
linearizing the spectrum for the spinous, we have 

ImG?(ki? — q,Lu — Q) = —n5{uj — Q + Ex + avpqz). 

Since, O ~ w, and q > q* for the overdamped mode, the 
constraint from the (5-function is always satisfied. After 
the angular integral, we get 



ImEf(fcF,w) 



f 



2B 



\ Trpo / Jo 



dnn 



dq- 



9* 



47r2f72 ' 



where momenta and frequencies are dimensionless in 
units of kp and aD, respectively. For fi ^ ui > E*, 
the leading contribution of the g-integral comes from 
q r-^ (ri)^/^ and therefore the infrared cutoff q* can be 
set to zero. But ior fl ^ lu < E*, q* > il^/^, and the 
lower cut-off at q* comes into play. We finally get. 



lmY.^{kF,Lu) 






yGnpoaD 



kp 



Gn'^poaDE* 



\lu\ > E* 



oj\ \uj\ < E* 



(40) 



Thus we find that above the infrared cutoff scale E* , the 
Kondo breakdown scenario, in which the conduction elec- 
trons interact with the critical hybridization fluctuations, 
provides a microscopic mechanism to obtain a marginal 
Fermi liquid in three dimensions. 

Next we evaluate the temperature dependence of the 
imaginary part of the self-energy at w = and on the 
Fermi surface. Denoting this as Iml]c(r), we get from 
Eq. 



ImEf (T) = 



q 



dn 



1 



, TT sinh(r2/r) 

X ImGf (ki. - q, -17). 



Imi9^(q,f]) 



(41) 



The evaluation of the above expression is very similar to 
the finite frequency case, except for T > E*, the thermal 
factors nsix) + np{x) = l/sinh(x) gives an additional 
logarithm which is cut off by E* . We get. 



ImEf (T) « - 



l3 

tup 



SirpoaD 

1.3 

fSjp 

6poaDE* 



Tlnl^ 
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T> E*, 



rpZ 



T <E* 



(42) 
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FIG. 5: The imaginary part of the conduction electron self- 
energy at tj = versus T from a numerical evaluation of 
Eq. (|4ip using Eq. (|25p . Note the approximate linear T be- 
havior, with a crossover to T^ behavior at very low T. 



In Fig. [HI we show a plot of this quantity from a nu- 
merical evaluation of Eq. (|4T|l using Eq. ([25l) . One can 
see the approximate linear T behavior except at very low 
temperatures, where one crosses over to a T^ behavior. 

Next we evaluate the temperature dependence of the 
resistivity 5p{T) = p{T) — /o(0). In order to proceed, 
we first need to address whether the transport life- 
time Ttr can be identified with the quasiparticle lifetime 
Tc(uj,T) ex [lmSc(tt',T)]^^, whose frequency and tem- 
perature dependences are given by Eqs. (gOl) and ([1^ . 
For this, it is useful to compare the Kondo-Heisenberg 
model with a single band model. In the latter case, the 
two lifetimes have a different temperature dependence be- 
cause the leading contribution to the self-energy comes 
from forward scattering processes with momentum trans- 
fer g ~ 0, but which are not effective in relaxing the cur- 
rent. As such, when vertex corrections are taken into 
account, r^" acquires an additional temperature depen- 
dence proportional to q^ ~ T^/^. However, this is not 
the case for the Kondo-Heisenberg model which has two 
bands, one of light conduction electrons and the other 
of heavy spinous. Due to the constraint of half filling 
[Eq. dH)], the spinon current operator Jfi = at every 
site i. Therefore, it is guaranteed by gauge invariance 
that a vertex correction involving the exchange of a sin- 
gle a boson [Fig. [SJa] , which involves an external spinon 
current operator, is identically zero^i. The first non-zero 
vertex correction involves the exchange of two a bosons, 
and we expect such a correction to be small by a factor 
of a. This can be understood as well in a Boltzmann ap- 
proach'^-^ where the transport vertex correction 1 — cos(0) 
gets replaced by 1 — acos(0), which is essentially unity 
since a ^ 1. 

Consequently, in the present theory, the transport life- 
time is proportional to the quasiparticle lifetime. The 
physical picture that emerges from the above discussion 
is that, when scattered from a a boson (c ^ / -I- cr), a 
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conduction electron transmutes into a spinon and relaxes 
its current in the bath of the spinous. More formally, the 
expression for the conductivity (ac) obtained from the 
current-current correlator in the Kubo formalism is given 

by 



-^(|)i:/;j[^-.,(^)][,.cf(.,. 



We write 



[Gf(k,w)] ' = c^ - efc - ReEH 



2T,iLu,T)-- 



(43) 



(44) 



where for {cJ,T) > E* 



[rci^,T)]- 



-1+^ 



2k p 



\3TTpoaD 



Tin 



2T 

'e* 



M 
2 



(45) 

Here r is an elastic scattering lifetime of the conduc- 
tion electrons due to impurities, and sets the scale of 
the temperature independent part of CTc- We linearize 
the dispersion of the conduction electrons and replace 
the momentum sum by an energy integral, and we fi- 
nally obtain ac{T) — [poVpTc{0,T)]/3. This implies for 
E* <T < aD, 



' 2T 
5p{T)^T\n[ — 



(46) 



Therefore, the scenario of the breakdown of the Kondo ef- 
fect captures one of the most enigmatic features of heavy 
fermion systems close to quantum criticality, namely the 
quasi-linear temperature dependence of the resisitivity 
observed for most compounds over a large range of tem- 
perature. For T < E* , the usual Fermi liquid result is 
recovered and Sp{T) oc T^. It is interesting to note that, 
the recovery of the Fermi liquid T^ behavior of resistivity 
below a finite temperature scale in the quantum critical 
regime of YbRh2Si2 has recently been reported^. Fi- 
nally, we note that for the same reason that equates the 
single particle and transport lifetimes, the electrical and 
thermal transport lifetimes are the same. 



VII. CONCLUSION 

To summarize, we studied the Kondo-Heisenberg 
model in three dimensions using a fermionic represen- 
tation for the localized spins. The mean-field phase dia- 
gram in the T — Jk plane, where Jk is the Kondo cou- 
pling, exhibits a quantum critical point that separates 
a uniform spin liquid phase from a heavy Fermi liquid 
phase. In the uniform spin liquid phase, the Kondo hy- 
bridization between the conduction band and the band of 
fermionic spinous that constitute the local moments van- 
ishes, thereby indicating that in this phase, the Kondo 
effect fails to occur. For a Kondo coupling larger than the 
critical value {Jk > JKa)i ^ heavy Fermi liquid ground 





(b) 



FIG. 6: Diagrams for the current-current correlator for the 
evaluation of the conductivity in the Kubo formalism. Solid 
dots indicate current vertices, (a) Involves conduction elec- 
trons with self-energy corrections [see Fig. |4] . This contribu- 
tion identifies the transport lifetime with the conduction elec- 
tron lifetime, (b) Vertex correction involving the exchange of 
one a boson. It vanishes because the spinon current is zero 
due to the local constraint [see Eq. 0]. 



state is established with finite hybridization between the 
bands. This implies that at the quantum critical point 
{Jk = Jk^), the lattice Kondo energy scale Tk vanishes, 
indicating the breakdown of the Kondo effect for cou- 
plings smaller than the finite value Jk^ ■ 

In general, the size of the (hot) "Fermi surface" of the 
spinon band is different from that of the conduction elec- 
trons, and we characterized their mismatch by a wavevec- 
tor q* = kpQ—kp, where kpo is the spinon Fermi wavevec- 
tor and kp is the conduction Fermi wavevector. As a con- 
sequence of this mismatch, we found that two mean field 
solutions are possible in the Fermi liquid phase. First, 
one with a uniform hybridization, which is stabilized 
when the two band masses have the same sign. This is the 
standard Kondo phase which appears in the mean-field 
pseudofermion description of the Kondo lattice. Second, 
a novel Kondo phase with the hybridization modulated 
in space with wavevector go ~ l-2<?*, which appears when 
the two band masses have opposite signs (i.e., one band is 
electron- like and the other hole- like). Conceptually, this 
phase is analogous to the LOFF state of superconduc- 
tivity, and is characterized by nodes in space where Tk 
is zero. In this paper, we did not examine the physical 
consequences of the modulated hybridization, which will 
be the topic for future work. For the uniform case, we 
showed that at the quantum critical point the single ion 
Kondo scale T^ is approximately equal to the Heiscn- 
berg coupling Jh- This demonstrates that the Kondo 
breakdown is a consequence of the competition between 
the Kondo energy scale and the magnetic energy scale, 
even though there is no long range magnetic order in the 
present formulation. 

Then, we studied the effect of the critical hybridization 
fluctuations (excitations of the order parameter a) associ- 
ated with the vanishing energy scale Tk on the thermody- 
namic and transport properties of the system. We found 
that, due to the mismatch q* , the critical fluctuations are 
affected by energy scales E* ~ [aD / {2Tr)] {q* /kp)^ and 
Ex ~ avpq* , where a '^ Jh/D is the ratio of the spinon 



bandwidth Jh and the conduction bandwidth D. The 
propagator for the critical modes has several asymptotic 
structures in different regimes of frequency-momentum 
space, out of which the following two are important and 
readily understood, (i) For momentum q < q*, the spec- 
trum of the critical fluctuations lies outside the interband 
particle-hole continuum and therefore their dynamics is 
undamped and is characterized by a dynamical exponent 
z = 2. (ii) For most of momentum space {q > q*), the 
spectrum of the critical modes lies within the particle- 
hole continuum, and therefore have overdamped dynam- 
ics with exponent z — 3 (Landau damping). The lead- 
ing contribution to thermodynamics and transport is al- 
most always governed by the latter asymptotic structure, 
in contrast to most Ginzburg-Landau approaches, where 
only the critical modes within l/^(r) of the ordering vec- 
tor {q — 0) are important"^*. Above the temperature scale 
E* ^ 1 mK, this gives rise to anomalous metallic be- 
haviour, such as a specific heat coefRcient that diverges 
logarithmically with temperature, and the inverse life- 
time of the conduction electrons which has a T In T tem- 
perature dependence. The latter is a consequence of the 
conduction electrons scattering with the critical bosons 
with the dynamical exponent z = 3, which in three 
dimensions provide a microscopic mechanism to obtain 
marginal Fermi liquid behaviour. Since the spinous do 
not carry current, but are effective in relaxing the current 
carried by the conduction electrons, the T-dependence of 
the inverse particle lifetime also gives rise to a T In T be- 
haviour of the resistivity. From a scaling point of view, 
in this regime the frequency SI of the critical fluctuations 
scale as J7 ~ g'^, where q is their momentum. For T < E*, 
however, the infrared cutoff q* prohibits the z — 3 scal- 
ing, and the leading T dependence of the specific heat 
coefficient and the resistivity are Fermi liquid like. 

The Kondo breakdown scenario is promising in that 
it can explain one of the least understood features of 
the heavy fermions near quantum criticality, namely the 
quasi-linear temperature dependence of the resistivity, 
and the existence of multiple energy scales, over decades 
of temperature. 
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one band is electron-like and the other hole-like, the mean 
field theory admits a solution where the Kondo hybridiza- 
tion is modulated in space. This is a consequence of the 
mismatch between the two Fermi surfaces, and concep- 
tually is analogous to the LOFF state of superconduc- 
tivity2£i2i. In the following, we choose the conduction 
band to be electron-like and the spinon band hole-like, 
and linearize their dispersions. This gives e^ — e for the 
dispersion of the conduction band, where e = VFik — kp), 
and e^ — —a{e — vpq*) for the dispersion of the spinon 
band. For this case, we evaluate the static interband po- 
larization nyc(q, 0), whose general expression is given by 
Eq. (jl9[) . Approximating the momentum summation by 



E-f 



D 1-1 

de / dz, 
-D J-i 



where the conduction bandwidth D enters as an ultravi- 
olet cutoff for the energy integral, we get. 



n/c(q,0) 



Po 



■In 



In 



avp |(g*)^ -' 



(Al) 



It is easy to see that the maximum of — n/c(q, 0) is at a 
finite wavevector go where 



go « l-2g* 



(A2) 



Therefore for Jk > Jk^ , where the critical value of the 
Kondo coupling is given by 



Jk^ 



n/,(go,0)=0. 



the Kondo boson condenses in the Fermi liquid phase 
at a finite wavevector go (i.e., (cgo) ^ 0). This implies 
that the Kondo hybridization is modulated, with nodes 
in space where Tk vanishes^^. 
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APPENDIX A 



2. Calculation of the mean field free energy 

In this part, we give the technical details for the eval- 
uation of the mean field free energy at T = 0. This can 
be written as 



Fmf 



k,i— a,6 



(-4) + ^ + # + ^, (A3) 
Jk Jh ^ 



1. Spatially modulated mean field solution 

Here we demonstrate that when the conduction elec- 
tron and spinon masses have opposite signs, i.e., when 



a.b 



where 9{x) is the Heaviside step function, and e^.' are 
given by Eq. ^. We replace J2k ~* Po J de, and from 
the solution of the equations e^' =0, we get 
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]s.,i=a,b 



Po 



PO 

2 



de 



-D 



de 



(1 + a) (e + 
(1 + a) (e- 



e 



1/2- 



where s = CTo/a^/^ and si = CTg/(2D) + ©(l/D^). We 
expand the expression under the square root in powers 
of (1/-D) and keep terms up to 0(1/D^), since higher 
orders contribute to O(o'o) ^^^ beyond which we neglect. 
Performing the e-integral, to 0(0-0) accuracy we get 



Fmf __ PoD 

N 



2 r ^2 



2 
1 — a \a 



2poJh 3 
1 



Poc^o 



PqJk 

PoG^ (1 - 4a + c 
a2i:)2 (i_^^) 



where a constant part has been ignored. Since a <C 1 , in 
the terms proportional to CTq ^'^d "'O' '^^ retain only the 
dominant a-dependence, and get Eq. ([11 



APPENDIX B 

1. Asymptotic structure of the Kondo boson 

In this appendix, we determine the leading frequency 
and momentum dependences of the propogator for the 
Kondo boson Dcr(q, i^n) in the quantum critical regime 
using Eqs. ^, ^^ and ^. Its leading frequency 
dependence is given by the first terms in Eqs. (|23ap - 
(|23ep . while the next term is determined by comparing 
the static {q/kpY term in Eq. (P0|) with the sub-leading 
terms of Eqs. (|23ap - (|23ep . The asymptotic structure of 
D„{(\, iQ,n) in different regimes of frequency and momen- 
tum are as follows: 

(1) VLn<^E* = [{aD)/{2T:)]{q*/kFf. In this frequency 
interval there are three sub-regimes depending on the 
magnitude of the momentum q. We get 
(a) q < qni = [^n/E,j,]kF (where E^ = avpq*), 



D-\q,iar,) 



-Po 



Ex 



l + -(l + a) 



Zi i^ 
Ex 



(b) qni <Cg<(7*, 

-D<;^(q,if^„) «Po 

(c) g* < q < kp, 

D^^{<l,i^n) ~ Po 



1(JL 

4 \kF 



4 [kp 



Zi i^ 
Ex 



TT \Qr: 



2 avpq 



(Bib) 



(Blc) 



e 
D 



{l-a)^(e+^-] +4^2 



(1 



1/2- 



^) +K 



r 



(2) E* <^^n '^ Ex- In this frequency interval there are 
four sub-regimes given by 
iB.)q^qn2 = i^n/Exy/^q*, 



D^'{<i,inn) W -po ( V^ 



1 / N i^n 

1+2(1 + ")^ 



(B2a) 



(b) qo2 < g < q* 



D„^{c[,iiln) w -po 



Ex 



1 f Q 



3 \q 



, (B2b) 



(c) g* « q < qna = kF[{q*nn)/{akpD)]^/\ 

\ avpq J 

(d) go3 < 9 < fci=^, 



z-sgn(f7„) H 

z (7 



L»^^(q,zn„) « Po 



4 I fcj. 



TT |fi„ 



2 aw^g 



(B2c) 



(B2d) 



(3) Ex '^ ^n ^ ctE). In this frequency range there are 
five sub-regimes given by 

(a) q^qni = q*iEx/nny/^, 



D-^q^in^) 



Po 



(l-a) 



(b) gn4 < g < g* , 

.-l/„ .0 ^ ., Po 



D-'{<i,in^) 



(l-a) 



In 



In 



Zi£^^ \ H/x 



H/j. J ll Lfi 



, (B3a) 



—irt„\ 1/9 



£^x / 6 \q 



(c) g* < q <C <7s:25 = kp[n,J{aD)] 
Po 



(Bla) i:'<,'(q,iaO 



In ^ +1 + *:t 

\ avpq / 2 



(l-a) 
(d) gns < g < %26 = fci.[f^„/(a£')]i/2^ 



(B3b) 

sgn(il„) 
(B3"'c) 



D~^{q,i^„) ^ Po 
(e) go6 < g < fcF, 



2 avpq 2 



1 <L 

4 I A:j. 



1, ,02 

(1 + a)- 



{avpqy 



TT 117,, 



2 avpq 



(B3d) 
(B3e) 
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2. Spectral response of the Kondo boson 

In the paper, several simplified expressions were used 
for the spectral response of the Kondo boson. Here, we 
give a more complete account. The fc polarization bub- 
ble has some similarities to the Lindhard function'^^ , but 
also differs from it in important respects. In particular, 
the particle-hole continuum of the Lindhard function ex- 
ists for all momenta, while this is not the case for the fc 
polarization as a result of the mismatch between the con- 
duction and spinon Fermi surfaces. Wc have performed 
numerical calculations including the full quadratic dis- 
persion of the fermions, but they are very similar to re- 
sults we present here that are based on Eq. (|22aP plus 
the static curvature correction (last term in Eq. (|20p ). 
The advantage of using Eq. (|22ap is that it is valid for 
arbitrarily small a. All results here are for the retarded 
response function at T=0. We confine our discussion to 
the case where both conduction and spinon bands have 
the same sign for the mass. 

We begin with the fc bubble at g = 



Renfcio,n) 



Po 



1 



In 



\n 



avpq 



a 



\^ 



vpq 



(B4) 



This expression contains two logarithmic singularities at 
the energies Ex = avpq* and vpq*, where q* = kpo ~ kp 
is the mismatch vector between the conduction and 
spinon Fermi surfaces. The imaginary part of Hfc is sim- 
ply a step function with value Trpo/{l — a) between these 
two energies. Note that Hfc is not symmetric around zero 
energy (only the sum of it with Ucf would be). We have 
chosen kpo > kp. For the reverse case, the singularities 
would flip to the other side of the frequency axis. The log 
singularity at E^ plays an important role. It guaranties 
that the Kondo boson propagator, D = Jx/(l-t- J^H/c), 
always has a pole between zero and E^ ■ This pole is un- 
damped since Imllfc is zero below E^. 

The general structure of Hfc can be appreciated from 
Fig. [3 where the various domains for the imaginary part 
are shown. Note that the imaginary part vanishes in the 
regime we label as z = 2. For the positive frequency side, 
this is a triangle in {q, V,) space bounded by (0, E^) and 
(q*, 0). For low frequencies appropriate for the dispersive 
peaks of IniD, it will be sufficient to expand Eq. (|22aP 
for small H.. When we do this, we find 



ReD-^ =6- 



Po^ Jj^.l? 



q 



2avpq \q-q*\ 



paq 



(B5) 



where S is the deviation from the quantum critical point 
(QCP) and the last term is the static curvature correc- 
tion. Below the kinematic boundary, ImD~^ is zero, so 
the zeros of Eq. (jBSp in this regime give the collective 
mode dispersion, which for (5 = is 



n,M/Ex^0.5iq*/kpy{q/qy/ln 



\q-q*\ 



(B6) 



We compare this in Fig.[S]to the expression where the log 
in Eq. (jBSp is expanded for small q/q*, the latter being 
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1 
q/q* 
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FIG. 7: Regimes of Hfc- The various dashed lines are the 
kinematic lines corresponding to the zeros of the arguments 
of the logs in Eq. (|22ap . The regime denoted z = 2 has 
Imllfc = whereas the regime z = 3 has Imllfc oc Q. 
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■ Eq. B6 

■Eq. 24 

Eq. B8 



Eq. 25 



z=2 




k/' z=3 



0.5 1 1.5 

q/q* 



FIG. 8: (Color online) Dispersion of ImD for 5 — and 
q* /kp = 0.1. The undamped {z = 2) dispersion is to the left 
of the kinematic boundary, Q/Ex = 1 — q/q* (marked by the 
nearly vertical line), whereas the damped (z = 3) response is 
to the right. The results based on Eq. (|22a|l closely follow the 
expressions of Eqs. (|B6|I and (|B8|I . The simpler Eqs. pi)l and 
(|25p are used in the analytic calculations, and are quite good 
except for q near q* . 



Eq. ([M|) . Note that formally, Eq. (|B6p vanishes as q goes 
to q*, but this is of no concern, since the mode intersects 
the kinematic boundary before this occurs, and thus it 
terminates at a finite energy, corresponding to c ~ 0.1 in 
Eq. (EH). 

Above the kinematic boundary, ImD~^ is non-zero. 
For q > q*, it is 



ImD- 



2avpq 



(B7) 



This leads to a pseudo-Lorentzian behavior for ImD. 
The location of the maximum of ImD, denoted as F, can 
be found upon differentiation with respect to fl, leading 
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FIG. 9: (Color online) Plots of ImD for positive (left) and 
negative (right) Q.. The quantum critical point {5 = 0) is 
shown on the top, away from this {5 = 1) is shown on the 
bottom. The z = 2 dispersion is not visible on the scale of this 
plot. Note the approximate (anti)syinmetry of the damped 
[z — 3) response at the QCP as compared to away. This 
damped dispersion at the QCP closely follows the analytic 
expression of Eq. (|B8|l . The intensity scale for the bottom 
plots are a factor of ten smaller than the top ones. 



highly asymmetric for non-zero S. As S increases, the 
most intense part of the damped response moves up the 
kinematic boundary fl/Ex = 1 — q/q* and approaches the 
log singularity at q = 0,n — E^- In Fig. [TOl the disper- 
sion of the ImD maxima is plotted for various S. Note the 
reversed magnon-like dispersion of the undamped modes 
and the approximate hnear q behavior of the damped 
modes for non-zero S. 
2 




FIG. 10: (Color online) Dispersion of the ImD maxima for 
S ranging from zero (bottom curve) to 1 (top curve). The 
undamped modes are to the left of the kinematic boundary 
(dashed line), the damped modes to the right. Note the re- 
versed magnon-like dispersion of the undamped modes and 
the approximate linear q behavior of the damped modes for 
non-zero 5. 



to 



T/E, = 0.5iq*/kF)^iq/q*)yJ7r^+lir 



2 \q + q*\ 



(B8) 



which is also plotted in Fig. [S] If instead, we ignore the 
n term in Eq. (jB5|) . we get Eq. (^5]) instead. The latter 
is a true Lorentzian, and its dispersion is plotted as well 
in Fig. [5] Although formally Eq. (|B8|) vanishes as q goes 
to q* , the actual results based on Eq. (|22ap do not, and 
we again find c ~ 0.1 in Eq. (PS)1 . 

We finish this discussion by showing in Fig. [5] ImD 
based on Eq. (|22ap for both positive and negative fl for 
two cases, the quantum critical point [6 — 0) and some- 
what away {6 = 1). The collective mode is not visible on 
the scale of this plot, but we note that it is only present on 
the positive frequency side. The damped response is ap- 
proximately (anti)symmetric in il for ^ = but becomes 



3. Free energy 

Here we compute the free energy due to the excita- 
tions of the Kondo boson, whose expression is given by 
Eq. ([27|) . and take into account all the different asymp- 
totic structures of the propagator Z?(j(q, i^n) which are 
given in Eqs. (jBl[) - (jB3p . The goal of this exercise is to 
prove that for all temperatures T < aD, the leading con- 
tribution comes from that part of the phase space where 
the boson is overdamped with dynamical exponent z = 3, 
and whose propagator is given by Eq. (j25p . 

(1) T < E*. Since for the leading T dependence we 
expect ri ~ T, in this temperature regime D^(q, i^n) has 
three asymptotic forms which are given in Eq. (|Bip . Ac- 
cordingly, we split the g-integral into three parts, namely 
q < qni, qni < q < q*, and q* < q < kp, and denote their 
contributions as Fia, Fu, and Fic respectively. Keeping 
only the leading terms for each sub-regime, we get 
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We note that, since T < E* , the leading temperature 
dependence is due to the z ~ 3 mode whose contribution 
is given by Eq. (|B9cp . and thus F « Fic- 
(2) E* < T < Ex- In this temperature regime, 
Z?CT(q, ifln) has four asymptotic forms which are given in 



Eq. (jB2p . Now we spht the g- integral into four parts, 
namely q < qn2, qm < q < q* , q* < q < qm and 
Qns < q < kp, and denote their contributions as -Pja, 
F2b, F2c and F2d, respectively. Once again, keeping only 
the leading terms for each sub-regime, we get 
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After comparing the various contributions above, once 
again we find that the leading temperature dependence 
is due to the z = 3 mode, whose contribution is given by 
Eq. (|B10d|) . and we have F w Fz^. 

(3) Ex < T < aD. In this temperature regime 
-DCT(q, iJ7„) has five asymptotic forms which are given in 



Eq. (jB3|l . Now we split the q-integral into five parts, 
namely q < go4, qn4 < q < q* , q* < q < qn5, 
qn5 < q < qne and q^^ < q < kp, and denote their 
contributions as E^ai ^36 j -Pscj ^3d and E^^ respectively. 
Once again, keeping only the leading terms for each sub- 
regime, we get 
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As before, we find that the leading temperature depen- 
dence is given by the z = 3 mode, whose contribution is 



given by Eq. (jBlle|) . and we have F ^ F-. 
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additional contribution to the conductivity via the loffe- 
Larkin composition rules^®. This contribution, though, is 
of order a, and therefore does not change our findings, 
since a <C 1—. 

^^ J. M. Ziman, Electrons and Phonons (Oxford Univ. Pr., 
London, 1960), p. 376-377. 

^■^ G. Knebel, R. Boursier, E. Hassinger, G. Lapertot, P. G. 
Niklowitz, A. Pourret, B. Salce, J. P. Sanchez, L Sheikin, 
P. Bonville, H. Harima, and J. Flouquet, J. Phys. Soc. 
Japan 75, 114709 (2006). 

^^ Note that the z — 2 modes make sub-leading contributions 
to all quantities but the crossover lines. They also make 
no contribution to the conduction electron self-energy and 
transport due to kinematic constraints. 

^^ For spherical Fermi surfaces, as many ordering wave vec- 
tors as allowed by lattice symmetry will condense, each 
with a modulus of go ■ 

^^ see e.g., A. L. Fetter, and J. D. Walecka, Quantum, Theory 
of Many-Particle Systems, (Dover Publications Inc., New 
York, 2003), p. 158-163. 



